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abstract 

As is well known, coordinates of D-branes are described by x matrices. From generic non- 
commuting matrices, it is difficult to extract physics, for example, the shape of the distribution 
of positions of D-branes. To overcome this problem, we generalize and elaborate on a simple 
prescription, first introduced by Hotta, Nishimura and Tsuchiya, which determines the most 
appropriate gauge to make the separation between diagonal components (D-brane positions) and 
off-diagonal components. This prescription makes it possible to extract the distribution of D- 
branes directly from matrices. We verify the power of it by applying it to Monte-Carlo simulations 
for various lower dimensional Yang-Mills matrix models. In particular, we detect the topology 
change of the D-brane bound state for a phase transition of a matrix model; the existence of 
this phase transition is expected from the gauge/gravity duality, and the pattern of the topology 
change is strikingly similar to the counterpart in the gravity side, the black hole/black string 
transition. We also propose a criterion, based on the behavior of the off-diagonal components, 
which determines when our prescription gives a sensible definition of D-brane positions. We 
provide numerical evidence that our criterion is satisfied for the typical distance between D- 
branes. For a supersymmetric model, positions of D-branes can be defined even at a shorter 
distance scale. The behavior of off-diagonal elements found in this analysis gives some support 
for previous studies of D-brane bound states. 
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1 Introduction 



The most characteristic property of D-branes is that each coordinate of N D-branes is promoted 
to an X matrix X^, rather than a set of numbers [1]. One usually identifies the diagonal 
components as the positions of D-branes in the conventional sense, and the A^(A^ — 1) off-diagonal 
components as degrees of freedom of open strings connecting them. However, these A^^ degrees 
of freedom mix under the U{N) gauge symmetry of the system of D-branes, and therefore there 
is inherent indeterminacy in this identification of the positions of D-branes. 

However, of course, one should not totally abandon the concept of positions for D-branes. For 
example, when D-branes are "far apart" , there should be a gauge where distances between A^ points 
specified by diagonal elements are large, and then off-diagonal elements should be dynamically 
suppressed by the mass-term induced from the potential term, 

- ^r[X^X^2 = Y^^xt, - X^.flXr/ + 0((off diagonal)^). (1) 

fiy^u i<j 

Although the U{N) symmetry allows us to make arbitrary unitary transformations to the almost 
diagonal matrices described above, it is clear that it makes much more sense to call diagonal 
elements in the gauge where off-diagonal elements are very small as positions of D-branes, rather 
than those in the general gauge. 

In more general situations it becomes more obscure whether a similar concept of positions 
of D-branes exists. However, for many applications, it would be useful if one could define the 
concept of positions for D-branes (in terms of A^ numbers) for general situations, since it is much 
easier to grasp intuitively. One application we have in mind is the gauge/gravity correspondence 
[21 [3]: the black hole geometry on the gravity side should be related to the shape of the D-brane 
bound state on the Yang-Mills side, and in order to study the latter it will be essential to define 
the positions of D-branes in an appropriate way. 

In this paper we use a simple prescription, introduced in [5] for a slightly different purpose, 
to define positions of D-branes. The essential idea is to maximize the diagonal elements, in an 
appropriate measure, by an unitary transformation. For example, for the zero-dimensional matrix 
model (the bosonic analogue of an effective action for D(— 1) branes) 

S = -^Tr[X„X,f, (2) 

where d-i- = li ■ " " ^D) are N x N Hermitian matrices, we perform the U{N) transformation 
such that the quantity 

D N 
/j=l i=l 

takes the maximum value. It is obvious that diagonal components in this gauge coincide with the 
D-brane positions discussed above for the (almost) mutually diagonalizable case. 

We analyze, using the above prescription (and its generalizations) for determination of D-brane 
positions, the results of Monte-Carlo simulations for various zero- and one-dimensional Yang-Mills 
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matrix modelqf|. Our aim is twofold: the first is to verify the power of this prescription. In 
particular, we will see that the method allows us to detect the topology change of the distribution 
of D-branes; this is consistent with the expectation from the gauge/gravity correspondence, and 
corresponds to the transition between black strings and black holes. |§ The second is to give a 
criterion to determine when the use of this prescription gives rise to a sensible definition of positions 
of D-branes, in other words, to give at least a partial answer to the important question, "when 
are the D-brane positions well-defined?" The basic idea is to see whether off-diagonal elements are 
governed by simple Gaussian distributions which are determined by diagonal elements; in such 
a case one can regard dynamics of off-diagonal elements as sub-dominant, compared to that of 
diagonal elements, justifying the separation between diagonal and off-diagonal elements. We find 
that this criterion is satisfied with a relatively short length scale, compared to the size of the whole 
D-brane distribution, which validates the use of our method to examine the shape of the D-brane 
bound state. 

This paper is organized as follows. In § [2] we introduce the matrix models we consider, and 
review what the gauge/gravity correspondence implies for qualitative behaviors of the models. 
These models are simplified versions of maximally supersymmetric Yang-Mills theories for D- 
branes; we also explain the relevance of the simplified models for our purpose. In ^ 13. II we explain 
our main ideas using the zero-dimensional matrix model as an example: we explain in detail the 
maximal diagonalization procedure introduced above and our criterion to determine when this 
procedure provides a sensible definition of D-brane positions. We also show the result of the 
maximal diagonalization for the standard example, the fuzzy sphere. The results of Monte-Carlo 
simulations for the diagonal and the off-diagonal elements are respectively given in § 13.21 and § 13.31 
In § m we study the bosonic matrix quantum mechanics. We first clarify its phase structure in 
§ 14. 1[ In § 14.21 we make an appropriate generalization of the maximal diagonalization to non- 
zero spacetime dimensions by combining it with the T-duality transformations for the Yang-Mills 
theory [5]. In § 14.31 we apply the method to the Monte-Carlo data: we read off the topology 
change in the shape of the D-brane bound state, a counterpart of the black hole/black string 
transition. In §[5] we show preliminary results for a supersymmetric matrix model; in this case we 
find that our criterion is satisfied even at shorter distances between D-branes. We also discuss the 
implication for the structure of the ground state. § E] is devoted to conclusion and discussions. 

2 Dual gravity description and black hole/black string transition 

In this section we introduce a class of matrix models in lower dimensional spacetime, which 
will be considered in this paper. We also review the gauge/gravity correspondence [21 [3] which 
will give us some expectations about qualitative behaviors of these models, in particular phase 
transitions corresponding to the transition between the black hole (BH) and the black string (BS) 
on the gravity side [H [71 [8l [9l [ini [HI [E] . The presentation here basically follows [7] . For review 
articles, see [13]. A world-sheet approach is proposed in [TJ], to understand this duality. 

^ Yang-Mills theories are effective theories for the low-energy dynamics of D-branes. For high energy processes 
one should also incorporate the higher excitation modes of open strings. Ideas discussed in this paper would be also 
useful even if these higher modes are included. 

® The phase transition itself was identified and studied previously, but it has not been shown that this phase 
transition is accompanied by the change in the shape (including information in transverse directions) of the D-brane 
bound state. 
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The models we consider are simplified versions of the models for Dp-branes at finite tempera- 
ture. Low energy dynamics of Dp-branes is described by supersymmetric Yang-Mills (SYM) in 
(1 -|- p)-dimension with the U{N) gauge group, which can be obtained from the ten-dimensional 
M = 1 SYM through the dimensional reduction [1] . By Euclideanizing the time coordinate and by 
compactifying it with a period (3, we obtain the finite temperature theory with the temperature 

Let us start with the maximally supersymmetric U{N) SYM theory in (1 -|- 1) dimension at 
finite temperature, 

S2dSYM = Jxj^ i^^u + 2(i?^X,)2 - [Xi, Xjf + (fermion)) , (4) 

where the temporal direction t and spatial directions x are compactified with length 1/Th and 
L, respectively. Our notations are as follows: Th is the temperature (the subscript H means it 
is identified with the Hawking temperature of the dual black brane), A is the 't Hooft coupling, 
Xj (I = !,••• ,8) are adjoint scalars. For fermionic fields the anti-periodic boundary condition 
should be taken along the temporal direction. 

In the high temperature limit of SYM, Th — > oo, only zero-modes in the temporal direction 
survive; in particular, all fermionic fields decouple. Then we obtain the bosonic matrix quantum 
mechanics 



N ['^ /I 1 

Sidfios = dx Tr i^-{D,X'pf - -[X'p, X'q 



(5) 



where P = (t, I), X'j = Xj, and X^. is the adjoint scalar field coming from the t-component of the 
gauge field. We will study models of this type in § HI We remark that it is far from clear whether 
one can use the gauge/gravity duality for this high temperature limit; for the high temperature 
limit of a D-brane bound state, effects of the higher excitation modes of open strings are not 
negligible. The reason for studying these high-temperature Yang-Mills (YM) theories such as ([5]) 
is that they are very useful to test our ideas because (i) Monte-Carlo simulations are much more 
tractable for them, and (ii) one can gain insight into their properties, in particular their phase 
structures, from corresponding low temperature theories (where one can use the gauge/gravity 
correspondence). These models are also interesting in their own rights. 

In the low temperature region 0, it is conjectured in [21 [3] that SYM theory (jj]) is dual to type 
IIB superstring theory on a certain background. The background is obtained by taking the near- 
horizon limit of the near-extremal black 1-brane solution in type IIB supergravity compactified 
on a circle. The solution in the string frame is \15\ [3] 



[/6 




(6) 



^ The dimensionless quantity X/T^ can be considered as an effective coupling constant, and it is customary to 
call the low temperature region also as the strong coupling region. 
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where di = 2^tt^, and < x < L parametrizes the direction. Thermodynamic quantities 
should be calculated in the Einstein frame defined by C^^"**^*" = e~''^/^G^*""^. The Hawking 
temperature is 

T. = A, (7) 



4:TTy/Xdl 

and the ADM energy and the entropy are respectively given by 

1 



E 
S 



27 . 3^5^2 

1 



24 . 37r5/2A3/2 



(8) 
(9) 



In order for the supergravity description to be valid, the string excitation and winding modes 
associated to the x-direction should be much heavier than the KK modes along S^, that is, 



1 



(10) 



Note that string loop corrections are negligible, because the string coupling e*^ is vanishingly small 
in the planar limit, where is taken to be large while A and Th are finite. 

By applying the T-duality transformation along the x-direction, we obtain the black string 
solutioqj, which is also called the smeared L'O-brane solution, in type IIA supergravity. 



ds' 



a 




dU' 



IF 



+ dx^ 



N 



mV/' 1 

) L 



(12) 



where x varies from to (27r)^/-L. Again the string coupling e'^' is vanishingly small in the planar 
limit. Thermodynamic quantities like the Hawking temperature, the ADM energy and the entropy 
are unchanged by the T-duality transformation. This type HA supergravity description is valid 
in a different parameter region; the winding modes are now negligible when 



The a' corrections are negligible when 



Tu«{. 



Th « 



(13) 



(14) 



which is the same condition as before. 

In this solution Z)0-branes are smeared and form BS winding on the compact direction. The 
essence of the BS/BH transition is to compare the free energy of BS with that of a BH solution 
with the same charge and temperature. In the BH solution, DO-branes clump to a small region. 



® In this paper, we call this solution as the black string. We do not call the solution to the IIB supergravity (|6} 
as the black string. 
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Such a solution in compact space can be approximated very well by the black 0-brane solution to 
IIA supergravity in noncompact space. (See [16], for example.) The solution is fTbl [3] 




2 



[/7/2 (l _ ^ 



N [U^ 



5/2 



(15) 



where do = 2^7r^/^r(7/2). This solution is conjectured to be dual to a (1 + 0)-dimensional SYM, 
and A, Th are the 't Hooft coupling and the temperature for the dual theory. [f| The Hawking 
temperature is unchanged by the T-duality transformation and A is related to that of (1 + 1)- 
dimensional SYM by 

. _ ^(1+1) 

-^(1+0) - — —■ (J^oj 
The action of the (0+l)-dimensional SYM is 

M C^/Th /II \ 

SidSYM = jj^ Tr i^-{DtXjf - -[Xj,Xj]^ + (fermion) j . (17) 

We will study an analogous model in § [5l In the high temperature regime the action (jl7p reduces 
to the bosonic matrix model ([2]) with D = 10. We study D = 3 and D = 4 cases in § [3l 

The BH/BS transition point can be identified by a careful comparison of the free energy of 
the two solutions. The critical temperature is 

^ (18) 



below which BS breaks down to BH. Here A is the 't Hooft coupling constant in (1 + l)-dimension. 
Note that not type JIB but type IIA supergravity is valid at Tc, comparing with equations (fT3]) and 
(|14p . validating our description of the transition using IIA supergravity. If we fix the temperature 
Th and vary L, then BS appears at large L. This is not counterintuitive, because if L is large 
then the compactification radius in the T-dualized picture is small; it is natural that BH ceases 
to exist if the compactification radius becomes smaller - eventually it winds on the compactified 
dimension and forms BS. 

The gauge theory dual of this transition was previously studied using the Wilson line operator 
W = P exp Jq dxAx^ , winding on the compactified direction, as the order parameter. The basic 
idea is as follows. As W is unitary, its eigenvalues lay on a unit circle on the complex plane. They 
correspond to the coordinates of L'O-branes along the compactified direction. In the BS phase, the 
eigenvalues should be distributed uniformly on the unit circle and hence the trace of the Wilson 
line should be zero. On the other hand, in the BH phase the eigenvalues are clumped in a small 



^ Recently this duality has been confirmed very precisely by studying the strong couphng regime directly using 
the Monte-Carlo simulation |17l 1181 [191 120j . The expectation value of the Wilson loop [19] and the internal energy 
[20j are consistent with their counterparts on the gravity side including a'-corrections. An interesting approach for 
understanding the duality has been proposed in [21j . 
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region on the unit circle and hence the trace should be non-zero. Thus, the transition from the 
BS phase to the BH phase can be identified with the breakdown of the global U{1) symmetry, 
Ax Ax + const., which implies W e^^W. On the gauge theory side, the breakdown of the 
^7(1) symmetry has not yet been studied for the low temperature regime, where the gauge/gravity 
correspondence should be valid. In [7] the bosonic matrix quantum mechanics has been studied 
as the high-temperature limit of the (1 + l)-dimensional model, and the transition, precisely 
characterized by the breakdown of the U{1) symmetry, is found. See |22j for an earlier work. It 
is natural to consider this transition as the continuation of the transition in the low-temperature 
regime predicted by the gauge/gravity correspondence. In |23j this transition has been studied 
further in detail and it turned out that there are actually two successive phase transitions. They 
can be characterized as "uniform black string to non-uniform black string" and "non-uniform 
black string to black hole" transitions; here (non-)uniformity refers to that in the compactified 
direction. At low temperature, supergravity analysis suggests there is only one transition, from 
uniform black string to black hole. It is an open problem to clarify the structure of the phase 
transitions at intermediate temperature. 

Monte-Carlo simulations for the maximally supersymmetric matrix models associated with 
D-branes would allow us to compare the results with the gravity dual. In this paper we shall only 
consider simplified versions of matrix models for D-branes, with fewer scalar and fermionic fields. 
The reason for this is that the simplified models are suitable for testing our ideas about defining 
positions of D-branes, because the simplified models are qualitatively similar to the full matrix 
models for D-branes, and Monte-Carlo simulations are more tractable for the simplified versionj^. 
For example, as we will see in § HI the bosonic matrix quantum mechanics with fewer scalar fields 
has the same phase structure as the one with 9 scalars [23], which is the high-temperature limit of 
the Dl-brane matrix model. Also, simulation results for 4-supercharge [261127) and 16-super charge 
[T71 [T8] matrix models are qualitatively the same. 



3 Matrix model in zero dimension and black hole geometry 

In this section we analyze the bosonic zero-dimensional matrix model ([2]). We first explain 
two main ideas in our paper: the maximal diagonalization procedure, first introduced by [1], and 
our criterion which determines when the concept of D-brane positions (in terms of N numbers, 

^" The technical obstacles for simulating the full matrix models for D-branes are as follows. First, Monte- 
Carlo simulations for the full model is much more computationally heavy. This problem is severe especially for 
supersymmetric models, because of the necessity of the calculation of the fermionic determinant. In the case of the 
bosonic models, the simulation cost itself does not depend so severely on the number of the matrices D. However, 
we expect that the convergence to the large- A'' distribution is slower at larger D, because the diagonal components 
are distributed in R^, and hence typically we have only N^^^ meshes for each direction. For example for DO- 
branes, say, we have D = 9, and if we want to divide each directions into four parts, then A'' should be 4^, which is 
hopelessly large. Second, supersymmetric matrix models with 16 supercharges have the notorious sign problem due 
to the pfafSan of the Dirac operator being complex. The sign problem potentially makes the construction of the 
distribution of the diagonal components very difficult. The reason is as follows. The notion of the distribution of the 
diagonal components fundamentally relies on the expectation that physics of the matrix model can be studied by 
considering typical samples. Usually typical samples are easily obtained by Monte-Carlo methods. However, if the 
sign problem is present, it can be difficult to obtain these typical samples, because there might be many spurious 
configurations whose effects simply cancel out. The sign problem does not exist for supersymmetric matrix models 
with four supercharges |25l [24] . 



7 



not N X N matrices) makes sense. We then analyze results of Monte-Carlo simulations according 
to these ideas. 



3.1 Basic ideas 

As explained in the previous section, at high temperature, supersymmetric matrix quantum 
mechanics describing DO-branes (|17|) reduces to the bosonic matrix model ^ with D = 10. This 
model with smaller values of D has been studied extensively in [3]. The action is 

N ^ 

S = -jYl Tr[X„X,f, (19) 

fl,U=l 

where /x, = 1, • • • ,D and X^^s are traceless Hermitian matrices. 

First, let us introduce the maximal diagonalization procedure [4]. According to the standard 
interpretation, diagonal elements of X^s correspond to positions of D-branes and off-diagonal 
elements correspond to degrees of freedom of open strings connecting two different D-branes [1]. 
But of course "diagonal" and "off-diagonal" are gauge-dependent notions - in order to obtain a 
unique definition of D-brane positions, we should choose an appropriate gauge. 

Let us consider the case where all branes are well separated. In this case, open strings are 
long and heavy, and hence not excited easily. If we neglect these very small excitations, we can 
simultaneously diagonalize the matrices describing the branes, and it is clear that one can identify 
the diagonal elements in this gauge as the positions of D-branes. Now let us gradually decrease 
the distances between the branes. Then we come to a regime where open string excitations are not 
negligible, but are still very small compared to typical values of diagonal elements. In this case, 
although in general gauges off-diagonal and diagonal elements are of the same order, in some gauges 
diagonal elements are much larger than the off-diagonal elements. It is still at least natural to call 
diagonal elements in these special gauges as positions of D-branes. This consideration leads us to 
maximally diagonalize matrices using an unitary matrix U = Umax G U (N) , where we maximize 
the following quantity. 



D N 

EE 

/^=l i=l 



ux„u^ 



(20) 



Clearly, if X^'s are simultaneously diagonalizable, ?7max^^f^max are diagonal. In general, J/max-'^^f^max 
are as close to simultaneously diagonalizable as possible. It is reasonable to define positions of 

D-branes as diagonal elements in this gauge, (UmaxXf^Umaxj ■ Note this notion of positions is 

V /a 

gauge-independent, as gauge-equivalent sets of matrices end up with the same diagonal compo- 
nents after the maximal diagonalization, if we neglect the unlikely case where there are several 
minima for the quantity (j20p . 

We wish to note that the measure for the maximal diagonality we employ, (|20p . although nat- 
ural, is certainly not the unique choice. One might also try to maximize Yl^=i X]j=i | 
for example. An advantage of our measure (j20p is that it respects the SO{D) rotational symme- 
try of the model. Another important property of (|20p is that maximizing the diagonal elements 
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amounts at the same time to minimizing the off-diagonal elements, because the matrix norm, 



D N 



D N „ D N 



fj.=l i=l fi=l i^j 



(21) 



is invariant under unitary transformations. We will revisit this issue in § 14.21 

The important next step is to consider when this maximal diagonalization procedure gives a 
sensible definition of positions of D-branes. As is clear from the above discussion, it is the behavior 
of the off-diagonal elements that is crucial for the consideration of when D-brane positions admit 
a sensible definition. Again, when all branes are well separated, certainly positions of D-branes 
should be well-defined. We wish to stress, however, that off-diagonal elements are non-zero, 
although very small, even in this case. To be precise, let us remind that the action can be written 
as 

^ = ^T.T.^^ti - X^jflXr/ + 0((off diagonal)^). (22) 

fi^u i<j 

When distances between branes X^^ — X^- are large, the first term suppresses the excitation of 
the off-diagonal elements: because of the Boltzmann factor e~'^, the distribution of off-diagonal 
elements should become a Gaussian with a very narrow width. The second term representing the 
higher order interactions is negligible because of this narrow width. 

The criterion we propose in this paper is motivated by this observation. We require first 
that the distribution of off-diagonal elements have an almost Gaussian form, and then compare 
the width of the Gaussian with that calculated from (|22p . We will say that the positions of D- 
branes are well-defined when the observed width of the Gaussian agrees well with the theoretical 
prediction from (j22p . Put in another way, we require that the higher-order interaction terms 
to be negligible and the off-diagonal elements is simply governed by the quadratic part of S. 
The reasoning behind this is that, when this criterion is satisfied, the dynamics of the off-diagonal 
elements is trivial in the sense that it can be inferred from the information of the diagonal elements 
without any difficulty. This justifies to make a separation, which is essential in defining positions of 
D-branes, between diagonal elements, containing interesting dynamics, and off-diagonal elements, 
playing a sub-dominant role in dynamics of the system. 

A more intuitive way of understanding our criterion is to think it as the condition when we 
can use the simple "D-brane-|-open string" picture (neglecting the open string interactions). If 
0((off diagonal)^) terms are negligible, transverse modes 0, which are orthogonal to Xa — Xjj, 
behave as the harmonic oscillator with mass ^fN\Xii — Xjj\. The (z, j)-components can naturally 
be identified with open strings stretching between two D-branes sitting at Xu and Xjj, as the 
mass of each mode is proportional to length of the corresponding open string, as it should be. 
If this relation does not hold, this would mean that the distance between the branes is so short, 
many strings are excited, and branes and strings interact so strongly that the relative positions 
of the branes cannot be determined precisely. 



As is clear from (|22p . we should only consider transverse modes: oflt-diagonal elements in the longitudinal 
direction do not have mass terms, as they are gauge degrees of freedom. In other words, there are no longitudinal 
oscillations of open strings. Furthermore, one can show that the condition of the maximal diagonality implies that 
the off-diagonal components in the longitudinal direction are always zero. 
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Figure 1: Distribution of "diagonal components" of fuzzy sphere. Spin j = 250. 

We stress that our criterion, which we beheve to be reasonably well-motivated, is not the 
only possible one, and there can be other criteria which are also useful. We also want to note 
highly dynamical nature of the problem for lower dimensional matrix models: in these models the 
vacuum expectation value of scalar fields cannot be fixed by putting it as boundary conditions. 
Instead, they should be dynamically generated, as emphasized also in |28| . 

As an illustration of the maximal diagonalization procedure, we apply the maximal diagonal- 
ization to a standard example, the fuzzy sphere. We consider three matrices 



where Jj is SU{2) generator of spin j = {N — l)/2. We have numerically determined the unitary 
matrix Umax which maximizes (pOj) . |ij We find that the distribution of the diagonal components 
after the maximal diagonalization approaches to a two-sphere of unit radius as j becomes large. 
(See Figured!) This is the expected result; the fuzzy sphere should be a spherical distribution of 
D{—1) branes. We note that typical matrices encountered in Monte-Carlo simulations are more 
non-commutative compared to the fuzzy sphere. One estimate for the non-commutativity is the 
ratio 



which is very small when N is large for the non-commutative sphere. For matrices we encounter 
in Monte-Carlo simulations, this ratio is typically of order one [11[29]. Other methods to extract 
shape from non-commutative matrices are proposed in [30] and |31| . 

In this case the result of the maximal diagonalization is not unique; because of the 5*0(3) rotational symmetry 
of the fuzzy sphere (which can be also realized by a part of unitary transformations) the maximum is degenerate. 
Of course, for generic configurations appearing in Monte-Carlo simulations there is no such ambiguity. 




(23) 




(24) 
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3.2 Distribution of the diagonal components — the shape of D-brane bound 
state 

Now we study the bosonic matrix model (|19|) . Let us start with the distribution of the diagonal 
components. The distribution is S'0(L')-symmetric as expected. The basic quantity is the radial 
density function where r = \Xii\, which is normalized so that / Attv"^ p{r)dr = 1 for D = 3 

and J 2TT^r^p{r)dr = 1 for D = 4. To determine p{r), we have collected many configurations of 
matrices (for example, 4000 samples for D = 3, N = 8) and evaluated the distribution. 
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Figure 2: Distribution of diagonal compo- 

, I -C? I / ^ ■ n o 11 Figure 3: r = iXul vs n(r) in D = 4 model, 

nents, r = |Ajj| vs p{r), m D = 3 model. ° ' ' '^^ ' 

For D = 3, an almost uniform rigid ball is formed, as we can see from Fig. [21 The unit of 
length is determined from the normalization of the action ()19p . Its radius is about 1.7. On the 
other hand, for D = A the distribution is more like a shell, a three-sphere with radius ~ 1.2, as 
is shown in Fig. [3l It is not clear whether or not the thickness of the shell goes to zero at large N. 
We note that since D-brane positions themselves have uncertainty as will be discussed in § 13.3^ 
even if the thickness becomes zero its physical implication is not clear. 



3.3 Distribution of the off-diagonal components and validity of the procedure 

In this subsection we study the behavior of off-diagonal components and confirm that our 
criterion discussed in ^ 13.11 holds for a wide range of distances between D-branes. 

First, we should make a few technical observations regarding the distribution of off-diagonal 
components. We start from observing that the maximally diagonal condition does not fix the gauge 
completely: the measure for the diagonality (j20|) is invariant under the C/(l)^ transformation, 
generated hy U = diag{e'^^^ , ■ ■ ■ ,e*^^). Under this U{1)^ transformation diagonal components 
Xii do not change while off-diagonal components Xij are multiplied by e'*-^*"^^^. Hence it is natural 
to study the distribution of the absolute value \Xij \ . The density function for off-diagonal elements 
then acquires the usual volume factor 27r|Xjj| in addition to the Gaussian weight. Because 

In our algorithm, we seek the maximum by moving around on the group space, where each step is a multipli- 
cation by a randomly generated SU (N) matrix. We have checked that two independent diagonalization procedures 
give the same result. The computational cost increases with A'^, and for A'^ = 64 and 96, this consistency check 
sometimes fails because of the limitation of resources; we expect the result is at least qualitatively correct. 
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of the overall factor in the action, it is also convenient to rescale off-diagonal components Xi. 



as XJf" 



NXij, so that y 



I vnew 



becomes 0(1). 



If we neglect the interaction terms in (j22p according to our criterion, the statistical distribution 



of y for fixed value of Xi, 



Xjj should behave as 



p{y) = 2a^yexp {-a^y^) , 



(25) 



with 



iXj, 



Xjj\ = d. 



(26) 



The overall constant is fixed by the normalization J dxp{x) = 1. We find that (j25p provides a 
good fit in a large region of \Xii — Xjj\, see Fig. [J] for example. In general the value of a differs 
from d but the ratio a/d as a function of d is close to 1 except for small d. This confirms our 
expectation that when D-branes are far apart, the distribution of off-diagonal elements is the 
Gaussian following from the first term in (|22p . At the short distance scale our criterion breaks 
down. As we will show in § [5l in the super symmetric case our criterion holds even for shorter 
distances. 



D=3, N=32 
D=3, N=16 
D-3, N=8 



D-3, N-32 
D-3, N=16 
D-3, N=8 



0.6 0.8 



Figure 5: The distribution of off-diagonal 
components p(y), with 1.00 < d < 1.05, 
hi D = 3 model. For N = 32, the ansatz 
(j25p with a = 1.51 provides a reasonable fit. 
Finite- corrections are rather large. 

As we can see from Fig. [6] and Fig. [TJ (j26p approximately holds when {Xa — Xjj\ > 1.5, which 
is comparable to the radius of the D-brane distribution. This fact strongly suggests that, at this 
scale, the relative positions of the D-branes can be determined precisely. In Fig. [8] and Fig. [9] 
we plot the distribution of distances between two points, P{d), normalized by the condition 
J dxP{x) = 1, which is consistent with the result of |24j . We can see that many pairs of branes are 
sufficiently separated so that their relative positions can be determined precisely. Thus, even in 
the bosonic model, resolution for the D-brane position is rather good and the D-brane distribution 

For each configuration we have measured distances for N{N — l)/2 pairs of points. The distances between two 
points belonging to different configurations are not measured. 



Figure 4: The distribution of off-diagonal 
components p{y), with 1.50 < d < 1.55, in 
D = 3 model. For N = 32, the ansatz ^ 
with a = 1.67 provides a reasonable fit. 
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shown in § 13.21 makes sense. This gives enough motivation to apply the maximal diagonalization 
procedure to study the BH/BS transition, which will be done in ^ 14.31 

As we can see by comparing Fig. H] and Fig. [SI finite- corrections are larger at the shorter 
distance scale. This is probably due to the finite- correction coming from the Faddeev-Popov 
determinant for the maximally diagonal gauge. We have not evaluated the Faddeev-Popov deter- 
minant in full generality. For the simple and important case where only one off-diagonal mode is 
excited, the Faddeev-Popov determinant can be evaluated and is proportional to ^ — -|- 
where d = {Xu — Xjj\, z = ^y\X^^\'^ + l^trl'^ ^-i^d Xtr stands for the relevant off-diagonal transverse 
mode of X"''^'^. For large A^, the first term dominates and one can simply neglect the determinant 
as it only gives a constant factor when considering the distribution of z for fixed d. For finite A^ 
and small d, z spreads to large values and hence the second and the third terms of the determinant 
should give large contributions. 

4 Matrix quantum mechanics and black hole/black string topol- 
ogy change 

The action of the bosonic matrix model in (0 -|- 1) dimension is 



Here D = dx — i[A, ■ ] is the gauge covariant derivative and Xi{i = !,••• , dg) are adjoint scalars. 
When dg = 9, this model is the high-temperature limit of (1 -|- l)-dimensional maximal super 
Yang-Mills theory. (We have absorbed factors such as XTjj in ([5]) by redefining fields and the 
coordinate x.) 

As is discussed in §[2l this model is related to the black hole/black string transition. We begin 
with the study of its phase structure in § 14.11 Next we generalize the maximal diagonalization 
procedure to higher dimensions in § 14.21 Then in § 14.31 we apply this procedure to the matrix 
quantum mechanics and see the topology change. We find that, at the transition [7] corresponding 
to the BH/BS transition [6], there is a topology change which strikingly resembles the counterpart 
on the gravity side |16| . 

4.1 Phase structure of matrix quantum mechanics 

In this section we study the phase structure of the bosonic matrix quantum mechanics. The 
argument is parallel to that in [23]; the only difference is in the number of scalars, dg. 

As is well-known, the model has a global U{1) symmetry, multiplying the Wilson loop TrW = 
■^TrP exp{i Jq dxAx), winding on the compactified direction, by a phase factor. This symmetry 
can be broken in the large- A^ limit. In order to detect the breakdown of this symmetry, we should 
measure (|TrVK|) rather than (TrW), because the overall phase factor is not fixed at finite- A^. 
As shown in [TJ [23] the global U{1) symmetry is broken when L is small. To study the detail 
of the transition, it is better to look at eigenvalues of W rather than TrW itself. Because W is 
unitary, its eigenvalues are written as e*^, where — vr < 9 < tt. We study the distribution of the 
phases p{9). Note that we fix the overall phase factor such that TrW becomes real. By looking 
at the distribution p{6), it was found that there are two successive transitions at L = Li and 




(27) 
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L = L2 < Li [23]. Above L = Li the phase distribution is uniform, p{9) = 1/{2tt), and hence the 
U{1) symmetry is not broken, (|TrH^|) = 0. Therefore, the theory is volume independent [321 133j 
in this region. Befow L = Li the distribution is not uniform. For L2 < L < Li, p{6) is not a 
constant and is non-zero everywhere. Befow L = L2 the distribution has a gap (i.e. p{9) becomes 
zero at a certain value of 6). It has been observed that the phase distribution below L = Li can 
be fitted [23] by an expression used for the Gross- Witten model [M], which is, 

— (l + -cos(9) (k>2) (28) 
27r \ K J 

for ungapped phase L2 < L < Li and 



- - sin2 - (k < 2) (29) 

for gapped phase L < L2. 

In the following we show the results for dg = 2 and dg = 3. In Figure [10] we show the phase 
distribution for dg = 2 when L is very close to L2. We find the gap emerges at L2 ~ (1.3)-^ ~ 0.8. 
This transition is of second order [23]. We note that the corresponding transition in the Gross- 
Witten model is of third order. 





Figure 10: Distribution of the Wilson line phase near transition point, dg = 2, N = 32, A = 12. 

In order to determine the value of Li, we utilize the L-dependence of k. For = 32 and 
A = 12, K"^ can be fitted by the straight line, k"^ ~ 3.31'^ - 3.7, at 1.2 < < 1.26, and 
we can estimate the value of Li, where k becomes infinity, to be around (1.12)-^ ~ 0.89. To 
determine the order of the transition, we fit the quantity J —jj:{Tr[Xi,Xj]'^)dx oc ^ log Z, where 
Z is the partition function, by an ansatz 

j -^(rr[A„ X,f)dx = a(L-i - 1.12)f + C, (30) 

where C ~ 0.94 is the value for L> Li. Then, for the data at 1.2 < L^^ < 1.26 for A^ = 32, A = 12, 
we perform a two parameter fit with respect to a and p, and we obtain p ~ 2.0, which means that 
the transition is of third order. 

For dg = 3, the behavior of the eigenvalues are qualitatively the same and, in this case, we 
find Li ~ (0.93)"^ ~ 1.08 and L2 ~ (1.1)"^ ^ 0.9. 
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4.2 Maximal diagonalization for matrix quantum mechanics: T-dual case 

As explained in § [2l we should go to the T-dualized picture in order to see the black hole /black 
string transition. In the matrix model, we should use Taylor's T-duality for Yang-Mills theories [5], 
which we shall briefly recall here. The starting point is the zero-dimensional bosonic YM, 

Sod = -\Tr[xP,X^]\ (31) 

where P,Q = (0, i) = (0; 1, • • • , ds). The next step is to divide into infinite number of blocks 
with the size A^, N^^^X^^; we obtain 

<^ = - —fr (X^ yQ - X^\ ( X^ X^ - x^ X^ ) C^?) 

4 \ mg qn -^mq-^qn) y^nr^rm ^nr^rm) ' \"'^/ 

where tr is a trace in each block and the summation over m,n, - ■ ■ is assumed. Then we impose 
the compactification constraint, 

^mn — ^(m-l)(n-l)' 

^mn = X(rn-l){n-l) '^)' 

^nn = X'^n-l){n-l) + ^^T-R • 1. (33) 

Under this constraint, the action (j32p should be divided by the (infinite) number of copies and 
re-expressed in terms of X^ = X^q. For n these matrices represent open strings between 
D-branes with the winding number n. For n = 0, diagonal components represent positions of 
D-branes and off-diagonal components represent open strings with no winding. By T-duality, the 
winding number is translated to the KK momentum: we choose variables after performing the 
T-dual transformation as 

A^ = Y^ e2™^^X°, y' = Yl e^™'^'^^;. (34) 

n n 

Then, the action becomes 

Sm = NrJJ dx tr Q(2?..y.)' - \[Yi,yjf^ ■ (35) 

Then, by rescaling 



we obtain 



Sid = N I dx tr i-{D^Y,f - -[Y„Y,f ) . (37) 







From these rules, we learn that open strings between different D-branes in the T-dualized 
picture correspond to the non-constant modes and the off-diagonal elements in the one-dimensional 
YM. Hence, it is natural to minimize non-constant modes and off-diagonal elements. We note 
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that we should distinguish here between the minimization of off-diagonal components (and non- 
constant modes) and the maximization of diagonal components, as one cannot use the conservation 
of the matrix norm by unitary transformations (|2ip . (We use the term "maximal diagonalization" 
throughout this paper although it is a slight abuse of the terminology. ) Actually, the minimization 
procedure gives sensible results, as we will see, whereas the maximization procedure turns out to 
give pathological results. One way to understand the difference is to note that the measure for 
the size of off-diagonal (and non-constant) elements and diagonal elements are positive. For the 
minimization procedure, this serves as a lower bound to stabilize the procedure, whereas the 
maximization procedure is more susceptible to pathology because of the absence of an upper 
bound. 

Another technical issue is that, to look at the distribution of branes, we have to take care 
of the effect of the translational symmetry, and fix the origin for branes. To fix the origin for a 
noncompact direction, we simply take X to be traceless. To fix the origin of the direction, 
we fix the global U{1) symmetry so that the Wilson loop rrPexp(z J dxAx{x)) becomes real and 
positive [23] . 

Our simulation is based on the non-lattice technique |26j . We generate configurations consist- 
ing of the gauge field in the static-diagonal gauge 

Ax = Y ■ diag{ai,- ■ ■ ,aN)-, L = i?"^/^, (38) 

and Fourier modes of adjoint scalars 

Yi{p) {p=-A,--- ,A). (39) 

Because the coordinate basis is more convenient for maximal diagonalization, we introduce 2A + 1 
lattice sites 

X = — L, ■■■ , 0, • • • , — L, (40) 

2A + 1 ' ' ' ' 2A + 1 ' ^ ^ 

and perform the Fourier transformation 

A 

nx) = ^ e'-'P'^/'^Mp). (41) 
p=-A 

The gauge field is mapped to the unitary link variable as 

U{x) = (iza(7(e^"i/(2A+i)^ . . . ^ ^ia^/{2A+i)^ ^42) 

After performing the unitary transformation in the coordinate basis to maximally diagonalize, we 
have to go back to the momentum basis to use the T-dualized picture. 
From the above considerations, it follows that one should minimize 



^-^Y.Tr[2-U{x)-U\x))+Y: 

X i 

The detailed derivation is given in Appendix [Al 



i,Re|x:(f/..(x)-l)| -p^^E^-c 



x)] 



17 



4.3 Black hole/black string topology change 



In this section we show the D-brane distribution, obtained by using the maximal diagonal- 
ization procedure introduced in the previous subsection. The off-diagonal components behaves 
similarly to the zero-dimensional model studied in § 13.31 and hence the distribution is meaningful. 
In Figure [TTl [T2] and [13] we plot diagonal components after the maximal diagonalization for the 
dg = 2 model. The x-axis represents the compactified direction and two edges ±1 should be 
identified. (In these plots we have rescaled the coordinates so that the period of the compactified 
direction becomes 2.) The y-axis represents the Xi direction. The X2 direction is projected out. 

When the radius R in the T-dual picture is small (i.e. L large) they form a uniform string, 
as we can see from Figure [TTJ If R becomes large, at some point the string begins to pinch; see 
Figure [T2j We may call it as the nonuniform string. Finally it is pinched off completely and 
becomes a squashed ball, as shown in Figure [T3j These shapes are strikingly similar to the black 
string and the black hole in gravity side [16]. As we increase R further, we have found that 
the distribution exhibits SO{ds + 1) symmetry, as expected. If one does not apply the maximal 
diagonalization and just plots diagonal components resulting from a Monte-Carlo simulation (in 
the static-diagonal gauge), these gradual change of the geometry is almost completely obscured. 
In particular, one cannot see the recovery of the spherical symmetry for large R. 

We have also studied the dg = 3 model, and found that the shapes of the bound states look 
similar. The difference is in their internal structures: for dg = 2 the bound state is uniformly 
filled {D^ X for the black string phase and for the black hole phase), and for dg = 3 the 
distribution is shell-like (S^ x for the black string phase and for the black hole phase). 




Figure 11: Distribution of di- 
agonal components for dg = 
2,N = 16, A = 8,L-i = 0.5. 
The "uniform black string" 
phase. 



Figure 12: Distribution of di- 
agonal components for dg = 
2,N = 16, A = 8,L~i = 
1.2. The "nonuniform black 
string" phase. 



Figure 13: Distribution of di- 
agonal components for dg = 
2, = 16, A = 8,L-i = 1.5. 
The "black hole" phase. 



5 Supersymmetric matrix models 

In this section we show preliminary simulation results on a supersymmetric matrix quantum 
mechanics. We adopt the non-lattice Euclidean path-integral method [26] . A Hamiltonian formu- 
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lation can be found in [35], and a mean field approximation is studied in [36]. As explained in 
§ m at low temperature the maximally symmetric model (|17p describes the black 0-brane in IIA 
supergravity. Here we study a simpler model, the four-supercharge analogue of PT|) : 

S = N 1^'^^ dt Tr (^^{DtXjf - ^[Xj,Xjf + ijDtij - ija'[Xj,i^]^ , (44) 

which is obtained from four-dimensional M = 1 SYM through the dimensional reduction. Here 
I, J runs through 1,2 and 3, are Pauli matrices, and tpa (« = 1)2) are N x N complex 
fermionic matrices. This model does not suffer from the sign problem and the simulation cost is 
less expensive, while we expect its behavior should at least qualitatively similar to that of (jl7p . 
To consider the finite temperature system we have to impose the anti-periodic boundary condition 
(a.p.b.c.) to the fermion. However, as pointed out in |17j . with the a.p.b.c. the system with small 
and finite T has an instability (eigenvalues of Xj run to infinity) , which makes the computation 
heavy as we need to take A^ to be reasonably large. This lead us to adopt the periodic boundary 
condition (p.b.c), as the p.b.c. and a.p.b.c should be the same in the T — > limit which we are 
interested in. 

Because the compactified direction is now temporal and the T-duality transformation is not 
performed, the maximal diagonalization procedure introduced in § 14.21 is not adequate. Rather, 
we maximize J dt (A^p) , that is, we maximally diagonalize scalars at each time slice, lj 

Numerically, at A^ = 6, the diagonal components we find are distributed like a ball. To determine 
the shape of the distribution reliably, it is also necessary to study larger A^ and see the A^- 
dependence. We would like to report on this issue in a future work. 

We should also introduce a slightly different criterion for determining the validity of the max- 
imal diagonalization. We again regard the diagonal and off-diagonal elements to be the position 
of the DO-branes and open strings, respectively. Our criterion is again to see when we can neglect 
higher order interaction terms (and accordingly can think of off-diagonal elements as dynamically 
insignificant). If interaction terms are negligible, and the diagonal elements are slowly- varying 
compared to the off-diagonal elements, each off-diagonal mode behaves as a harmonic oscillator. 
Quanta of this harmonic oscillator correspond to open strings and have energies proportional to 
their lengths. As we are interested in the low temperature limit, we furthermore require that the 
open string is not excited at all, the wave function being of the Gaussian form. This is the minimal 
non-commutativity; off diagonal components cannot vanish because of the zero point oscillation. 

Indeed, we find that the distribution of y = \/iV|Ajj| can be fitted by the probability distri- 
bution function for the zero-point oscillation 

p{y) = Aaye-^'^y' (45) 

in a wide region, and a is very close to d = \Xii{t) — Xjj{t)\. In Fig. [H] we plot a/d against d. 
In Fig. [T5l we plotted a/d for fixed d against T. At short distance, the data begins to depart 
from the ansatz (j45p . We suspect the reason to be corrections in the Faddeev-Popov determinant, 

A similar maxinial diagonalization (and the criterion for its validity), which does not count the gauge field (or 
the scalar field) in the temporal direction exists for for the zero-dimensional matrix model studied in § [S] This is 
more suitable than that given there (which treat the gauge field in the temporal direction in an equal footing as 
other scalar fields) if we regard this system as the high temperature limit of the DO-brane system. The resulting 
D-brane distributions from these different prescriptions do not differ qualitatively. 
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which are not neghgible for finite N as we have seen in § I3.3t to study d < 1.5 reUably, we have 
to take larger N. 

The apphcability of the ansatz (|45p is at least an indication of the vahdity of the assumption 
that the DO-branes (diagonal components) are more slowly-varying than the open strings (off- 
diagonal components). We note that this assumption is crucial for previous studies |37^ [38] for 
D-brane bound states, using various analytic and approximate methods. For another approach 
to bound states, see [Ml HOI HIl IMJ- For a discussion of bound states for the zero-dimensional 
matrix model, see 



N=G,A=4,T-1.5 I — 
N-6,A=e,T-0.8 ^- 
A=8. T-0.4 
N-e,A=ie,T-0.2 f 
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t A mi J.- I 1 ■ A OTTOAA X ■ Figure 15: Plot of aid versus T at d = 
l^igure 14: ihe ratio a/ a m 4 bUhY matrix ' 

, V. ■ ■ 1 xx J ■ x J 1.55,2.05 and 2.55 in 4 SUSY matrix quan- 

quantum mechanics is plotted against a. 

tum mechanics. 

From these plots, we can see that the notion of the D-brane distribution becomes better at 
smaller T. A more interesting possibility suggested by these plots is that a/d = 1 holds for 
any (i at T = in the planar limit. Let us elaborate on this possibility and the general tendency 
that our criterion holds better at lower temperature. As we are taking the T — > limit, it 
is expected that the off-diagonal elements are in ground states, not in excited states. What is 
surprising is that the effect of the higher-order interaction terms seems to be negligible even for 
small distances between D-branes, because for small distances the width of the wave function (for 
each off-diagonal element) is large and therefore one expects that the higher order terms in the 
action also contribute. The fact that this absence of the effect of interaction terms occurs for 
supersymmetric models and not for bosonic models suggests that this might be a consequence of 
the familiar cancellation between bosonic and fermionic degrees of freedom. It is a very interesting 
problem to understand the mechanism of this phenomenon. Of course, to obtain a rigid conclusion 
we also have to study the model more thoroughly for larger N and smaller T. We would like to 
report on these issues in a future publication. 



This behavior is in sharp contrast with that in the bosonic model shown in Fig. 1161 In the bosonic case, 
the temperature dependence disappears at low temperature. This is as expected, because for the bosonic model, 
the global U{1) symmetry (multiplying the Wilson loop by a phase factor) is recovered at low temperature, which 
implies that the large- A'^ dynamics does not depend on the volume (in this case the temperature) [32ll33j . as is well 
known. The consistency with this general result makes us rather confident that we are taking a sufficiently large 
N to neglect finite- A'' effects. In the supersymmetric case the global U{1) symmetry is recovered only at T = [44] 
and hence the dependence of T does not disappear. 
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Figure 16: The ratio a/d in the bosonic matrix quantum mechanics with 3 scalars is plotted 
against d. 

It might be worthwhile to compare this with an approach |45l I46| . centering around the 
Berenstein's conjecture [l5] on the strongly coupled = 4 SYM on M x S^. As this approach 
focuses on dynamics of the diagonal elements, neglecting that of the off-diagonal elements, it is 
closely related to the themes of this paper. In this approach one has an additional dimensionful 
parameter, the mass for scalar fields, originating from the curvature of the . Because of this, 
under a crucial assumption explained below, the strong coupling limit implies very large mass 
terms for the off-diagonal elements. This in turn implies that the width of the wave function (or 
the Gaussian) is very small, which automatically validates the neglection of interaction terms. 
This is similar to, but has a crucial difference from what seems to be happening in our model: 
for our models, the widths of the Gaussian distributions are generated dynamically, and are finite 
instead of vanishingly small. Nonetheless, also in our model higher order interaction terms seem 
not to have effect in dynamics in the low temperature regime. 

The assumption mentioned above, underlying the analysis of = 4 SYM on M x S^, is that 
the typical length scale for distances between D-branes should be determined basically by the 
additional dimensionful parameter (the size of S'^) and should not depend too much on the coupling 
constant, which is also dimensionful. In [29], a mechanism for explaining this behavior is proposed, 
which is also closely related to the behavior we find numerically for T — > 0. The mechanism is 
based on an instability of the D-brane bound state in finite-iV supersymmetric matrix models 
(without mass-terms for scalars), where distances between D-branes become arbitrarily large. Ell 
In the presence of the mass-term for scalars, there would be a competition between this instability 
and the attractive force from the mass-term. The balance between these competing effects would 
be achieved at the point where the coordinates of D-branes take values of the order of the length 
scale determined by the mass-term, thus validating the assumption. One can partially justify the 
existence of this instability by using the one-loop approximation |27l BH] which is valid once the 
distances between D-branes become large. Our numerical results suggest that one might be able 
to extend this result even for shorter distances. 

This instability, which is of interest by itself, is suppressed severely for large A*'; the bound state we have 
discussed in this section might be called as metastable in this sense. This instability is numerically found in [17] , 
and presumably is the same as that first discussed in [47] , 
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As also noted above, the low temperature behavior we find numerically, namely the absence of 
the efi'ect of higher order interaction terms for a rather short distance scale, suggests that physics 
in the low temperature regime might be understood by an one-loop approximation. We hope that 
our numerical results provide a guide to such an improved perturbation theory. 

6 Conclusion and Discussions 

In this paper we have proposed a simple procedure to determine the shape of the D-brane 
bound state from super Yang-Mills theories. The strategy is to maximally diagonalize the matrices 
representing the collective coordinates of the D-branes. We have introduced a criterion which 
determines when the notion of this D-brane position makes sense. We have tested the use of these 
ideas by Monte-Carlo simulations. As the maximally supersymmetric SYM are computationally 
hard, we have studied simpler models- zero-dimensional and one-dimensional bosonic YM and 
one-dimensional SYM with four supercharges. In the bosonic models, the criterion is satisfied for 
the typical distance scale between D-branes. The size of the whole distribution of D-branes is 
a few times larger than this length scale below which one cannot talk about D-brane positions. 
This is natural because a theory of quantum gravity would have a minimal length scale and very 
hot black holes, presumably described by the bosonic matrix models, would be objects of the 
size comparable to this scale. We also observed the topology change of the D-brane distribution 
corresponding to the BH/BS transition, by combining the above procedure and the T-duality 
prescription for YM. In the supersymmetric model, the notion of the D-brane position seems to 
be precise even at shorter distance scales. 

There are many future directions. Most straightforward one is to try to understand the 
gauge/gravity correspondence, and furthermore, to study the black hole physics via the corre- 
spondence. We have already seen the analogue of the BH/BS transition in the bosonic model. 
It is in principle straightforward to repeat the similar analysis for (l-l-l)-dimensional SYM. The 
lattice formulation for (l-l-l)-dimensional SYM has been developed recently [l9] and already 
Monte-Carlo simulations have been performed for four-supercharge SU{2) gauge theory [50]. By 
carrying out simulations for larger and larger supercharges the BH/BS transition can be di- 
rectly studied. At low but finite temperature, we expect to see even the stringy correction to the 
transition. As we have reviewed in § [2l the pattern of the phase transition changes at intermediate 
temperature - although there are two successive phase transitions at high temperature, the dual 
gravity analysis predicts there is only one transition at low temperature. It should be possible 
to see how high and low temperature regions are interpolated. At present it is still numerically 
challenging, but will be possible in near future. Also it will be interesting to study (1-1-2)- and 
(l-|-3)-dimensional SYM compactified on torus, which are expected to have richer phase diagrams 
|5H I52j. A generalization to curved spacetime is also interesting. For this purpose, the T-duality 
procedure on curved spacetime [53] and a construction [5l] of the four-dimensional = 4 planar 
SYM on from the BMN matrix model [55] might be useful. 

Another interesting direction is to apply our techniques to study noncommutative spaces and 
super membranes [56]. As we can see from Figured] the maximal diagonalization procedure can 
be a powerful tool to read off the shape of noncommutative manifolds. In this case, off-diagonal 
components are very small, although non-zero, and the shape is determined precisely. For example, 
we might be able to see the shape of the higher-genus membranes (fuzzy surfaces) constructed 
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along the line of ^T^. For the case of membranes, the embedding of the topology in the matrices 
has been discussed in [31], and it will be interesting to see the relation with the current approach. 
The study of the ground state of the matrix models for string/M theory |56l I58| . which 
presumably is associated with a different type of large N limit than the planar limit, is also an 
important direction. 

We wish to stress that the rather high applicability of the Gaussian distribution for the off- 
diagonal elements are, although natural, far from trivial. This Gaussian behavior suggests an 
existence of an improved perturbation theory (presumably some kind of mean-field approximations 
incorporating the information of the distribution of the diagonal elements) and its validity at the 
one-loop level. Usually one does not expect that the perturbation theory should be useful, as the 
't Hooft coupling constant A is dimensionful and we cannot compete it with the kinetic term for 
the mass-less case we are interested in. (For matrix models with a mass-term the ratio between 
the mass and the 't Hooft coupling constant gives a good estimate about the applicability of the 
perturbation scheme. 0) It is also interesting that the validity of the perturbation is associated 
with each pair of D-branes rather than with the whole configuration. As shown in § [5] for the 
supersymmetric matrix quantum mechanics at low temperature, the one-loop approximation is 
valid for very short distance scales, and the one-loop effective action provides a precise description. 
This observation may support a proposal in [21], that a critical exponent 2.8 of the energy density 
of the black hole, ~ 7AT'^-^, predicted by using the supergravity dual [3], might be understood 
from the one-loop effective action [60]. It will be nice if we can reproduce not only the exponent 
but also the overall factor 7.4. More importantly, the distribution of D-branes given by our method 
should correspond somehow to the metric of the black hole. Clarification of this correspondence 
would give us a good clue to understand the gauge/gravity duality directly. We hope that our 
method provides a useful mean to organize data and extract physics from it. 
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A Details of the maximal diagonalization in matrix quantum me- 
chanics in T-dual picture 

In this appendix we explain why we should minimize (|43p in the maximal diagonalization 
procedure. 

One might think that the dimensionful parameter T can be used to construct a dimension-less expansion 
parameter. Although this is useful for some systems, the expansion parameter is relevant only for the decoupling 
of Kaluza-KIein modes in the temporal direction, and does not validate the usual perturbation expansion based on 
Gaussian integrals for our matrix models without the mass terms. In [59], a hybrid approach, in which KK modes 
are integrated out perturbatively and zero modes are treated by Monte-Carlo simulations, is studied. 
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Because 



(j-x+a 
i J dx'A{ 



x') 



JaA{x) 



U{x) + U\x) c^2-a^A{xf, 



we have 



(46) 



(47) 



P ij 



Here a = L/(2A + 1) is the lattice spacing. Also, the zero- mode of diagonal components can be 
approximated as 



L 



L 



(48) 



Therefore, the contribution to the quantity we wish to minimize from the off-diagonal and nonzero 
modes can be written as 



P V 



2 ~ 2A + 1 



Y.Tr{2- U{x) - C/t (x)) + ^Y.^Ay. (Udx) - 1) | 

X i \ X ) 



Expressions for adjoint scalars are also similar. We have 

Y,TrY{xf = {2A + l)Y,Y.\Mp) 



P ij 



and 



y(o) = —^Yy{x), 
^ ^ 2A + 1 ^ ^ ^' 

X 

and hence the total contribution from adjoint scalar is 



(50) 



(51) 
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(52) 



As TrY{x) is gauge invariant, we need to evaluate only the second term. In summary, we should 
minimize 



2A + 1 



L2 



^ Tr (2 - U{x) - U^{x)] + 



i,Re|^(f/.(x)-l)| -j^j^{y^Ux)^ 



Note that first term does not appear if we perform the maximization of diagonal components 
instead of the minimization of off-diagonal components. 
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